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Abstract We seek for self-similar solutions describing the time-dependent evolu- 
tion of self-gravity systems with either spherical symmetry or axisymmetric disk 
geometry. By assuming self-similar variable x = r j at where a is isothermal sound 
speed we find self-similar solutions extending from the initial instant i = to 
the final stage i — > cxd using standard semi-analytical methods. Different types 
of solutions are constructed, which describe overall expansion or collapse, enve- 
lope expansion with core collapse (EECC), the formation of central rotationally 
supported quasi-equilibrium disk as well as shocks. Though infinitely many, these 
self-similarity solutions have similar asymptotic behaviors which may impose di- 
agnosis on the velocity and density structures in astrophysical systems. 
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Far away from initial and boundary conditions, dynamical evolution of fluid systems may lead 
to self-similar phases with physical variable profiles shape-invariant and magnitudes properly 
scaled (Sedov 1959; Landau & Lifshitz 1959). Similarity methods, which transform partial dif- 
ferential equations (PDEs) to ordinary differential equations (ODEs), greatly simplify nonlinear 
problems. For spherical systems, several well-known self-similar solutions were found to describe 
the collapse of an isothermal cloud in the context of star formation (Larson 1969; Pension 1969; 
Shu 1977; Hunter 1977; Whitworth & Summers 1985). The first is the LP solution independently 
found by Larson (1969) and Pensten (1969); the second is the "expansion wave collapse" solution 
(EWCS) found by Shu (1977), who also discovered other solutions with central point mass from 
the initial instant i = to the final stage t -l~oo; Hunter (1977) discovered infinitely many 
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discrete analytic solutions within the pre-catastrophic period which share with the LP solution 

that they arc regular at a; ^ 0~ although some of these solutions can not be continued in the 
post-catastrophic period from t = to t ^ +00. Hunter's results were extended by Whitworth 
& Summers (1985, hereafter WS), who allowed weak discontinuities across the sonic critical line 
and hence constructed continuous bands instead of Hunter's discrete solutions. But their solu- 
tions with weak discontinuities were criticized by Hunter (1986) that are only of mathematical 
interest and are physically unstable. Recently, Lou & Shen (2004, hereafter LS04) re-examined 
this classic problem and derived new solutions in the 'semi-complete space' (0 < i < -|-oo), in 
contrast to the 'complete space' (—00 <t< +00) taken by Hunter (1977) and WS. 

The realistic situation may deviate from the spherical collapse. In essence, the presence 
of rotation or magnetic fields will cause the collapsing system to be more and more flattened 
(e.g., Nakamura, Hanawa, & Nakano 1995, hereafter NHN) so that the system will later evolve 
disk-like. Models considering self-similar evolutions in a thin disk have been studied in many 
papers (e.g., NHN; Li & Shu 1997; Saigo & Hanawa 1998; Krasnopolsky & Konigl 2002) with or 
with out rotation and magnetic fields. Axisymmetry is assumed in the above papers to simplify 
the problem. 

In this work, wc plan to establish the whole solution structure for an isothermal fluid system 
with either spherical symmetry or axisymmetric disk geometry within the period from t = 
to t — > 00. In comparison with the previous solutions found by other authors, we particularly 
emphasize the existence of such solutions that the innermost part (core) is c;oIlapsing towards 
free-fall while simultaneously the outer part (envelope) is expanding and approaching constant 
wind. However, the overall collapse region is expanding in a self-similar manner. These envelope 
expansion with core collapse (EECC) solutions exist in the spherical and non-rotating disk cases 
but not the rotating disk case (see Section 3.2). Furthermore, as an non-trivial extension we can 
incorporate purely azimuthal magnetic fields in the disk cases (details are not discussed here). 



2 SPHERICAL CASES 



The basic formalism for the problem with spherical symmetry is well-established by many 
authors. By the definition of self-similar variable x = r/at where r, a, t are the radial coordi- 
nate, isothermal sound speed and time coordinate respectively, the reduced ordinary differential 
equations (ODEs) are as follows 
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where a(x), v{x), m{x) are the reduced density, radial bulk velocity and enclosed mass, respec- 
tively. The physical density p, radial bulk velocity u and enclosed mass M are obtained by the 
similarity transformation 
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m{x). 



(2) 



By expanding the solutions in power serials at large and small x; and with careful treatment 

of the sonic critical line one may construct various types of self-similar solutions which are 
present in LS04. In addition to the shock-free solutions present in LS04, we can also construct 
shocked similarity solutions, which are applied to the star- forming regions and H 11 regions 
(Shen & Lou 2004). 
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Among various solutions, we are particularly interested in such solutions, for which the 
innermost part (core) is collapsing and approaching free-fall towards the center while simulta- 
neously the outer part (envelope) is expanding and approaching a constant wind. The core-mass 
accretion rate is constant. The shock which connects the inner part and the outer part is prop- 
agating outward at a constant speed. 



3 AXISYMMETRIC DISK CASES 

For axisymmetric disk cases we introduce similarity transformation: 

r aaix) / \ ^ 

x = — IJ(r, t) = , u(r, i) = av{x) , Ca = 

at ZttGi 

M{r, t) = —m{x) , j(r, t) = Pa^tm{x) , $(r, t) = a^(/)(, 

where 



(3) 



— i r ' -^(^^ = / (l + X^-?Xcos^)V. ' W 

and a, w, q, m, cj) are reduced quantities which are dimensionless functions of similarity 
independent-variable x and S, u, Ca, M, j, $ are surface mass density, radial bulk velocity, 
Alfvenic speed, enclosed mass, specific angular momentum and gravitational potential respec- 
tively. The specific angular momentum is proportional to the enclosed mass because their ratio 
is assumed to be spatially uniform at the initial instant t = and conserved during the succes- 
sive evolution (Li & Shu 1997; Saigo & Hanawa 1998). Using these reduced variables we derive 
the ordinary differential equations (ODEs) as 



[{x-vr-l-q^]^^^f~^-^-^^ + p-^, (6) 

a ax x-^ X 2x X — V 

where / = d(j)/dx together with the following auxiliary equations: 

m = xa{x — v) , (8) 

and equation (4). Equations l@}-ljHI) are the full set of ODEs to be solved. In particular, equation 
©-(CI) are simultaneously solved for reduced density, reduced radial velocity and Alfvenic Mach 
number. In the monopole approximation, the reduced gravitational force is approximated as 
/ ~ m/x^. 

There exists one exact solution for the ODE set, which is the rotational equilibrium state 
of a (magnetized) singular isothermal disk (MSID): 

1 + D^-Qy2 1 + D^-Qy2 2 , 

V = , a ^ , / = , m — (1 + D ~ Q /2)x , 

X X j-g^ 

^= 1 + D^^-Qy2 ' g^constant^Q. 
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Without magnetic fields (g = 0) and in the monopole approximation, we reduce the basic 
equation set (Q-® into 



2 11^" a{x - v)[l - P"^ a{x - v)] - 1 

2 Til'^" a[l- (3^a{x-v)]-{x-v) 

[(X — V) — 1\ — = (x — V) , m = xa(x — v) , 

a ax X 

with the asymptotic solution when x — > +oo as 

l-A + p^A"^ ^(1-/32^2) 2F2 + (A-1-/?2a2)(A-4) 
X 2x'^ 6x^ ' 

A A(l- A + (3^A^) 
a ^ 1 —: , m Ax . 

X 2x-^ 

For non-rotating disks we simply set P — m the above equations. 

3.1 Non-rotating Disks Without Magnetic Fields 

For non-rotating P — cases, we have asymptotic solutions when a; ^ as either 



(10) 



(11) 



v ^ 

which stands for free-fall, or else 
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which stands for core-expansion. And the expansions of solutions near the sonic critical line 
X — V = 1 are 

v = x^-l + C^{x-x^) + ... , a = 1 + {C^ - l/x^){x - x^) + ... , (14) 

where 

^ 1 T yi + 2(l/xg - 1/x.) ^ ^^^^ 

with -minus and -plus signs denoting type 1 and type 2 solutions in reference to the spherical 
case (Hunter 1977; WS; LS04). Figure 1 shows several example solutions for the non-rotating 
case. 



3.2 Rotating Disks Without Magnetic Fields 

For rotating /3 > cases, notice from equation (5) that the term of centrifugal force diverges as 
x~^ if m mo, faster than other terms. So for rotating disks the material will not fall freely 
to the center as indicated by p2|l . Instead, we have asymptotic solutions when a; — > as cither 
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or the quasi-equilibrium asymptotic solution 
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Self-similar solutions 





disk-EWCS 



Fig. 1 Example solutions in non-rotating disks without magnetic fields. Left: Example 
solutions of v(x) versus similarity variable x. The dash-dotted line is the sonic critical 
line X — V = 1. The light solid line is one of the Shu-like solutions and the heavy solid 
line is the disk-EWCS. The dotted line is one of the disk-EECC solutions which do not 
cross the sonic critical line. The dashed line is the disk-LP solution. Right: Example 
solutions of a{x) versus similarity variable x. Notations and line-types are the same 
as left. 



where 

A^ -^"-f^ >-i, n^l±^^>u (18) 

and 77 is a constant parameter. Asymptotic solution (|17() is actually the quasi-equilibrium state 
with the leading order as the exact equilibrium SID solution and is rotating supersonically with 
rotational Mach number D exceeds unity. 

By the isothermal shock jump condition (Tsai & Hsu 1995; Saigo & Hanawa 1998; Shu et 
al. 2002) 

{Vd - Xs){Vu - Xs) ^ I , —^{vd~Xsf, (19) 

we can connect the inner quasi-equilibrium solution to the outer envelope which approaches 
constant winds or inflows or breezes. In the above equations Xs is the shock location fixed in 
the similarity coordinates which suggests that the shock moves outward at constant speed x^a, 
i.e., the physical radius of the shock discontinuity increases linearly with time. 



4 DISCUSSION 

The major difference of our self-similar solutions from those by other authors is that our so- 
lutions include not only those with constant inflow speeds or being quasi-static at large radii 
but also those with constant winds at large radii. Spatially, the novel solutions describe both 
collapsing part and expanding part, connected with or without shocks. Temporally, the fluid 
element at all radii is decelerated from expansion to collapse and the stagnation points as well 
as shock fronts propagate outward at constant speeds. The density and velocity profiles at 
large or small radii have well-defined power-law behaviors on radius, which thus can be used 
as useful input parameters in numerical radiation transfer codes. Another important aspect of 
these self-similar solutions is the core-mass accretion rate, which is constant. So the mass of 
the central object (either a proto-star or a black hole) increases linearly with time. We have 
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extended the range of the core- mass accretion rate, which is only of specific values for several 
known self-similar solutions (i.e., for the EWCS, the core- mass accretion rate mo = 0.975). 

The astrophysical applications of these self-similar solutions may appear in various circum- 
stances: systems involving accretions and outflows. To be more specific, the evolution of young 
stellar objects (YSOs), H ii regions around OB stars and planetary nebulae (PNe) may expe- 
rience certain stages which can be described by these self-similar solutions (Shen & Lou 2004). 
By incorporating more realistic ingredients such as the polytropic equation of state, the toroidal 
environment and even general relativity, these solutions originally derived in isothermal spheres 
can be generalized to describe more complicated cases (e.g., core-collapse supernovae). In ref- 
erence to our novel EECC solutions, we may propose the following scenario: the envelope at 
large radii is initially expanding due to external heating (i.e., from supernovae or evolved stars 
near star-forming regions); while the core region begins to collapse due to gravitational insta- 
bilities; the successive evolution of the core-envelope system may evolve into a self-similar state 
described by one of the EECC solutions. Such a situation may be applicable to the protostcUar 
evolution where core accretion and large-scale outflows may concur. For planetary nebulae, it is 
interesting to note that in the classic interacting stellar wind (ISW) model, a much faster stellar 
wind from the core catches up with a slower dense wind the remnant of the asymptotic giant 
branch phase-and together they form shocks; while our shocked EECC solutions imply that 
while the shocked expanding envelope is creating the observed planetary nebula, the core (i.e., 
the proto-whitc dwarf) is continuously accreting material. If it can accrete enough material to 
exceed the Chandrasekhar limit, there might be the possibility of igniting a Type la supernova 
explosion, without the requirement for a companion star! 
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